alpha-nucleus potentials for the neutron-deficient p nuclei 
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a-nucleus potentials are one important ingredient for the understanding of the nucleosynthesis 
of heavy neutron-deficient p nuclei in the astrophysical 7-process where these p nuclei are produced 
by a series of (7,11), (7,p), and (7,0) reactions. I present an improved a-nucleus potential at the 
astrophysically relevant sub-Coulomb energies which is derived from the analysis of a decay data 
and from a previously established systematic behavior of double-folding potentials. 
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I. INTRODUCTION 



The bulk of the heavy nuclei (A > 100) has been synthesized by neutron capture in the s- and r-process. However, 
most of the rare neutron-deficient nuclei with A > 100 cannot be produced by neutron capture. The main production 
mechanism for these so-called p nuclei is photodisintegration in the astrophysical 7-process by (7,n), (7,p), and (7,0) 
reactions of heavy seed nuclei from the s- and r-process. A list of the neutron-deficient p nuclei from 74 Se to 196 Hg 
^T^ \ can be found in Table 1 of Ref. . Typical parameters for the 7-process are temperatures of 2 < Tg < 3 (Tg is the 
temperature in GK), densities of about 10 6 g/cm 3 , and time scales in the order of one second. Several astrophysical 
sites for the 7-process have been proposed, and the oxygen- and neon-rich layers of type II supernovae seem to be 
a good candidate. However, there has been no definite conclusion reached yet with respect to the astrophysical site 
where the 7-process occurs. Details about the astrophysical scenarios can be found in the reviews by Lambert [Q, 
Arnould and Takahashi 0, Wallerstein et al. fj, and in Refs. 

Almost no experimental data exist for the cross sections of the 7-induced reactions at astrophysically relevant 
energies. Therefore, all reaction rates have been derived theoretically using statistical model calculations. One 
striking example is 146 Sm which is a potential chronometer for the 7-process. The production ratio of 146 Sm and 
144 Sm depends sensitively on the (7,11) and (7,a) cross sections at the branching nucleus 148 Gd, and it has been shown 
that especially the (7, a) cross section can be calculated only if a reliable a-nucleus potential is available. Predictions 
from different potentials differ by one order of magnitude |9|~fl2"||. The need for improved a-nucleus potentials for 
\ astrophysical calculations has been pointed out in p3jJT4] . 

Systematic a-nucleus potentials have been presented in several papers (e.g. Refs. [p^|-p7|), and recently these studies 
have been extended to astrophysically relevant energies (p]|l4|] . Usually potentials are derived from scattering data. 
• • ■ However, at the astrophysically relevant energies below the Coulomb barrier it is difficult to derive the potential from 
the experimental data unambiguously (see e.g. [p~Tj ) . Alternatively the cluster model provides another possibility 
to determine a-nucleus potentials by an adjustment of the a-nucleus potential to the bound state properties of the 
nucleus (A + 4) = A ® a p8|-p0[| . The half-lives of many a emitters have been calculated in [ pl]j22]| using a specially 
shaped potential, but this potential had to be modified to describe elastic scattering for 208 Pb j23||. It has been shown 
in Refs. that a systematic double- folding potential is able to reproduce both the bound state properties of 

212 Po = 20"Pb (g) a and elastic a scattering of 208 Pb. Finally, folding potentials give an excellent description of the 
experimental data on 144 Sm at 20MeV @, but again the potential which was used for the calculation of the a decay 



data 1 21 1 is not able to reproduce the precision scattering data. 

The basic idea of this work is to determine the a-nucleus potential at astrophysically relevant energies of about 
5 to 12MeV which is in the energy gap between the bound state potentials and the scattering potentials. I have 
chosen double folding potentials because of the small numbers of adjustable parameters. A systematic behavior of the 
strength of double folding potentials at higher energies is already given in Ref. [jlTj . In this work I will analyze bound 
state properties. First, I will briefly present the method in Sect. |lj then I will give results for the neutron-deficient 



potentials (Sect. IV) 



a emitting p nuclei in Sect. [II, and finally I will give an outlook on possible further improvements of the a-nucleus 
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II. FOLDING POTENTIALS AND a DECAY 



The a decay of nuclei is a clear proof that the wave function of the a emitting nucleus A + 4 has a non-negligible 
component A ® a, whereas low-energy elastic scattering is described using a pure A® a wave function. Therefore, an 
effective a-nucleus potential should describe simultaneously the half- life of the a emitter A + 4 and elastic a scattering 
of the nucleus A. In the astrophysically relevant energy region the a-nucleus potential is accessible mainly from the 
decay data because elastic scattering is dominated by the Coulomb interaction. For the analysis of the a decay data 
I apply the semi-classical model of Ref. ^5|, and the nuclear potentials V/v(r) are calculated by the double- folding 
procedure: 



V}v(r) = \V F (r) = A / / pp{r P ) p T (r T ) v eS (E, p = p P + p T ,s = \r + rp - rr\) d r P d r T 



(2.1) 



where pp, px are the densities of projectile resp. target, and u e ff is the effective nucleon-nucleon interaction taken in the 
well-established DD M3Y parametrization |26| , |27| . Details about the folding procedure can be found in Refs. 0,0, 
the folding integral (2.1) has been calculated using the code DFOLD | p8| . The strength of the folding potential is 
adjusted by the usual strength parameter A with A ~ 1.1 — 1.3 leading to volume integrals Jp per interacting nucleon 
pair of about 300 to 350MeVfm 3 . Jp is defined by 



Jp 



-in 



A P A 



PSlT Jo 



Vjv(r) r 2 dr 



(2.2) 



Note that in the discussion of volume integrals J usually the negative sign is neglected; also in this paper all J values 
are negative. 

The densities of the nuclei have been derived from the experimentally known charge density distributions |29| and 
assuming identical proton and neutron distributions. For nuclei where no experimental charge density distribution is 
available (i) the density distribution of the closest neighboring isotope was used with an adjusted radius parameter 
R ~ A 1 / 3 ( 186 0s, 182 W, 170 Yb, 150 Gd), (ii) the average between two neighboring stable isotopes was used ( 146 Sm) 
and (iii) the average of 138 Ba and 142 Nd was used for 140 Ce. 

The total potential is given by the sum of the nuclear potential Vn{t) and the Coulomb potential Vc(r): 



V(r)=V N (r) + V c (r) 



(2.3) 



The Coulomb potential is taken in the usual form of a homogeneously charged sphere where the Coulomb radius Rc 
has been chosen identically with the rms radius of the folding potential V F ■ 

The potential strength parameter A was adjusted to the energy of the a particle in the a emitter (A + 4) = A <g) a. 
The number of nodes of the bound state wave function was taken from the Wildcrmuth condition 



Q = 2N + L = ^(2n« + h) = ^< 



(2.4) 



where Q is the number of oscillator quanta, N is the number of nodes and L the relative angular momentum of the 
a-core wave function, and qi = 2n.j + li are the corresponding quantum numbers of the nuclcons in the a cluster. I 
have taken q = 4 for 50 < Z, N < 82, q = 5 for 82 < Z, N < 126, and q = 6 for N > 126 where Z and N are the 
proton and neutron number of the daughter nucleus (see also Table |). 
The a decay width T a is given by the following formulae pa| : 



T« = PF— exp 
4p 



-2 / k{r)dr 
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with the preformation factor P, the normalization factor F 
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(2.6) 



and the wave number k(r) 



k{r) = 



j \E-V(r)\ 



(2.7) 
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/x is the reduced mass and E is the decay energy of the a decay which was taken from the computer hies based on the 
mass table of Ref. |3Cj. The r» are the classical turning points. For + — > + s-wave decay the inner turning point is 
at n = 0. T2 varies from about 7 to 9fm, and varies from 45 up to about 90 fm. The decay width r Q is related to 
the half- life by the well-known relation T a = fi, In 2/^/2- 

It has to be pointed out that the preformation factor P should be smaller than unity because the simple two-body 
model assumes that the ground state wave function of the a emitter A + 4 contains a pure A® a configuration. The 
decay width in this model therefore always overestimates the experimental decay width. I determine the preformation 
factor P from the ratio between the calculated and the experimental half-lives |nj . A strong A® a cluster component 
is expected for nuclei A with magic proton and/or neutron numbers, and indeed the calculations show increased values 
for P around N = 126 ( 208 Pb) and N = 82 (e.g. 140 Ce) (see Table |). 

A preformation factor of P = 1 as used in |2lJ seems to be the consequence of the specially shaped cosh potential 
of that work. As an example I compare the potentials V(r) — Vm{t) + Vc(r) from this work and from for the 
system 190 Pt = 186 Os ® a in Fig.g. The rms radius of the potential from pl| is significantly smaller (r rms — 5.58 fm) 
than the rms radius of the folding potential (r rms = 5.97 fm). Therefore, the Coulomb barrier in |y| is significantly 
higher and the calculated half-lives in J2l| are roughly one order of magnitude larger than in this work. Note that in 
pl[ the potential was adjusted only to decay properties with the assumption P = 1 whereas in this work an effective 
potential is presented which is designed to describe decay properties and scattering wave functions and which leads 
to realistic preformation factors P. 



III. RESULTS FOR NEUTRON-DEFICIENT a EMITTERS 

The results of the calculations are summarized in Table | and shown in Fig. ||. One important result is that the 
strength parameters A and the volume integrals Jr for all a emitters show only small variations over the analyzed 
mass region A > 140. This means that the a-nucleus potential is well defined at very low energies. As expected from 
the systematic study in fl7j| , for the light system 8 Be = 4 He eg) a a much higher volume integral is required. 

The preformation factors P systematically increase to smaller masses with local maxima around the magic neutron 
numbers N = 82 and N = 126. The very high value of P = 65% for 8 Be is not surprising because of the well-known 
a cluster structure of this nucleus. 

One further exception is found for Hf = 170 Yb ® a with P = 62.8%. However, this surprisingly large value 
reduces to P = 3.13 % if the energy E = 2.584 MeV from plj is taken which was derived from the measured a energy 
and corrected for recoil and atomic effects instead of E = 2.4948 MeV from [p0| . On the other hand, the uncertainty 
of the measured half-life of 174 Hf is also 20 %, and a previous experiment gives a half-life which is roughly a factor 
of 2 higher |plj . From these calculations I have strong evidence that there is an inconsistency in the system 174 Hf = 
170 Yb <g> a between the measured half-life, the a energy, and the masses from the mass table |p0|| . 



IV. CONCLUSIONS AND OUTLOOK 




The real part of the optical potential is well defined by the systematic study of scattering data at higher energies 
above about 20 MeV and by the adjustment of the potential to the bound state properties at very low energies (see 
V . An interpolation between these energy regions leads to the recommended volume integral which is shown in 
as full line. A Gaussian parametrization is applied to Jr: 

J R (E) = J Rfi x exp [-(E - E ) 2 /A 2 } (4.1) 

with Jr q = 350 MeV fm 3 , Eq — 30 MeV, and A = 75 MeV. This interpolation is valid from very low energies up to 
about 40 MeV. 

The energy dependence of Jr at low energies is AJr/AE w 1.7MeVfm 3 /MeV, in agreement with AJr/AE = 
1 - 2MeVfm 3 /MeV and somewhat larger than AJr/AE = 0.71 Me Vfm 3 /MeV ftj. The uncertainty of an 
interpolated value of Jr is significantly smaller than lOMeVfm 3 corresponding to about 3 % in the interesting energy 
range around 10 MeV. 

Whereas the real part of the potential is well defined, no experimental information is available for the imaginary 
part of the potential at very low energies. However, it has been shown that transmission coefficients and cross sections 
in statistical model calculations depend sensitively on the volume integral J/ and even the shape of the imaginary 
part of the potential. Fig. || shows also the volume integrals Ji of the imaginary part of the potential for the same 
nuclei together with a parametrization from [ p2[ . A different parametrization of the energy dependence of J/ has 
been presented in Km. Because of the very similar behavior of 90 Zr and 144 Sm further information could be obtained 
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from a scattering in the A rs 100 region where the energy region between 10 and 20MeV might be accessible for 
high-precision scattering experiments. Note that in the 144 Sm case an analysis of scattering data below E = 20MeV 
Jll| is very difficult because of the dominating Coulomb interaction. 

New experiments in the A w 100 region are planned. If carried out with sufficient precision, the predictions of this 
work for the real part of the potential can be tested and new information on the imaginary part can be derived. 
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FIG. 1. Comparison of the potentials Vir) = Vzv(r) + Vc(r) from this work (folding potential, full line) and from ]2l|] 
(specially shaped cosh potential, dashed line) for the system 190 Pt = 186 Os ® a. The decay energy E is indicated by a dotted 
line. Note the significantly higher Coulomb barrier in fell] compared to the folding potential. 
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FIG. 2. The volume integrals Jr (upper diagram) and the preformation factors P (lower diagram) are shown for several 
neutron-deficient a emitters. The variation of the volume integrals Jr with A is small, whereas the preformation factor P 
increases to lower mass numbers A. Local maxima for P can be found around A = 144 and A — 208 corresponding to the 
magic neutron numbers N = 82 and N = 126. For 174 Hf = 170 Yb ® a the value derived from E = 2.584 MeV §l| is shown; 
the value derived from E = 2.4948 MeV @ exceeds the scale of the diagram and is indicated by an arrow (see text and Table 
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FIG. 3. Volume integrals Jr (upper diagram) derived from a decay data (this work) and from elastic scattering [jl7j|llj for 
the nuclei 90 Zr, 144 Sm, and 208 Pb. The full line shows the recommended interpolation for the astrophysically relevant energy 
region. The lower diagram shows the volume integral Ji of the imaginary part of the potential for the same nuclei together 
with Brown- Rho parametrizations j32] (dashed and dotted lines). This figure is combined from data of ]17[|ll|]. 



TABLE I. Experimental and theoretical decay widths of several neutron-deficient a emitters (parent). The light a emitter 
8 Be has been added. 



parent 


daughter 


E (MeV) 


Q 




A 


J fl (MeVfm 3 ) 


Ti? 2 c (y) 


7T/ X 2 P (y) 


P(%) 


!ua Po 


2U8p b 


8.9541 


22 


1 


241 


328.4 


9.20 x 10 _a s 


(2.99±0.02) x 10~ Y s 


3.08 a 


210p o 


206p b 


5.4075 


20 


1 


148 


304.5 


2.35 x 10' 3 


0.37886±0.00001 


0.62 


208p o 


204p b 


5215.5 


20 


1 


157 


306.8 


3.27 x 10" 2 


2.898±0.002 


1.13 


190 Pt 


186 0s 


3.2495 


18 


1 


067 


285.5 


3.01 x 10 10 


(6.5±0.3) x 10 11 


4.63 


186 0s 




2.8220 


18 


1 


078 


291.9 


9.08 x 10 13 


(2.0±1.1) x 10 15 


4.54 


174 Hf 


170 Yb 


2.4948 


18 


1 


116 


309.3 


1.26 x 10 15 


(2.0±0.4) x 10 15 


62.8 b 


154 Dy 


150 Gd 


2.9466 


18 


1 


125 


309.1 


1.16 x 10 5 


(3.0±1.5) x 10 6 


3.86 


152 Gd 


148 Sm 


2.2046 


18 


1 


146 


309.7 


1.42 x 10 13 


(1.08±0.08) x 10 14 


13.1 


150 Gd 


146 Sm 


2.8089 


18 


1 


141 


306.1 


1.73 x 10 5 


(1.79±0.08) x 10 6 


9.65 


148 Gd 


144 Sm 


3.2712 


18 


1 


159 


311.2 


6.20 


74.6±3.0 


8.31 


148 Sm 


144 Nd 


1.9860 


18 


1 


123 


303.2 


1.08 x 10 15 


(8±2) x 10 15 


13.5 


146 Sm 


142 Nd 


2.5289 


18 


1 


138 


307.1 


9.53 x 10 6 


(1.03±0.05) x 10 s 


9.25 


144 Nd 


140 Ce 


1.9052 


18 


1 


147 


307.9 


4.58 x 10 14 


(2.29±0.16) x 10 15 


20.0 


8 Be 


4 He 


0.0919 


4 


1 


624 


444.2 


4.34 x 10~ 17 s 


(6.71±1.68) x 10" 17 s 


64.6 



a Minor differences between the preformation factor P in Ref. |2j| and this work are due to the different choice of the Coulomb 

radius: R c = 1.2 fm x A^ 3 = 7.11 fm |m| and R c = 6.099 fm = r rms (this workV 

b A value of P = 3.13% is achieved with the energy E = 2.584 MeV from Ref. Q (see text). 
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